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Abstract 



The shell model embedded in the continuum is applied to calculate the astrophysical 
S-factor and the reaction rate for the radiative proton capture reaction ^^F(p,7)^^Ne. 
[ The dominant contribution to the cross-section at very low energies is due to Ml 

' transitions J[ = 2'^ —>■ Jf = 2f whose magnitude is controlled by a weakly bound 

• 22 state at the excitation energy Ex = 3.62 MeV. 

o ■ 

C3 ' Key words: Continuum shell model; Radiative capture reactions; Spectroscopic 

factors. 



1 Introduction 



Radiative capture reactions play a significant role in nuclear astrophysics being 
one of the important processes in stellar nucleosynthesis [1,2]. Nova explosions 
in binary star systems occurring due to mass transfer to a white dwarf star 
from its companion is an important site for the CNO cycle. A thermonuclear 
runway results from the accreted matter accumulated on the outer layers of 
white dwarf star, during which ^^F could be produced during hydrogen burning 
of the CNO elements. At this juncture the ^^F(p,7)^^Ne radiative capture 
reaction becomes very important in deciding the fate of ^^F and the subsequent 
routes of nucleosynthesis. The proton capture rate of ^^F is also important 
because in competition with the ^^F-/9-decay rate it can provide an alternate 
path from the hot CNO cycle to the rapid proton (rp) burning process via 
the "^''0(a,p)^'^F(p,7)^^Ne(/3"'', i^) ^^F(p,7)^^Ne reaction sequence. A smaller 
^^F(p,7)^^Ne reaction rate compared to the ^''F-/3-decay rate at characteristic 
nova temperature and densities would favour ^^O enrichment via the reaction 
sequence ^'^F{P^, z/)^'''0(p,a)"'^^N(p,7)-^^0, whose /?-decay to ""^^N could be an 
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explanation for the overabundance of ^^N in nova ejecta. On the other hand 
a large ^^F(p,7)^^Ne reaction rate would alter the ^^p^iTp abundance ratio 
due to the activation of the ^^F(p,7)^^Ne(/3"'", z/)^^F reaction chain. A good 
knowledge of the ^^F(p,7)^^Ne reaction rate is also thought to be important 
for understanding the astrophysical phenomenon of x-ray bursts. 

It was first predicted in Ref. [3] that the 3+, strong s-wave resonance in ^^Ne 
would occur at an excitation energy = 4.328 MeV with the width F = 
5 keV. This low energy 3"'" was expected to substantially contribute to the 
^^F(p,7)^^Ne reaction rate above stellar temperatures greater than 0.2 GK, 
which would be important for nova explosions. Since then there have been 
several attempts [4,5,6] to experimentally measure the exact position of the 
3"'" level in ^^Ne and it was only recently [7,8] that this has been unambiguously 
put at an excitation energy Er^ = 4.524 MeV (center of mass energy (cm.) 
Er = 600 keV) and width F = 18 keV. 

The direct capture part of the ^^F(p,7)^^Ne reaction has been studied largely 
in a potential model framework following mainly the work of Rolfs [9]. A 
more microscopic approach was the extended two-cluster model (ETCM) [10], 
using the generator coordinate method, which has also been applied to neutron 
rich Carbon isotopes [11]. One of the main advantages of this method is that 
the resonant and non-resonant part of the reaction need not be calculated 
separately. However the configuration space which is used to describe the 
many-body states is restricted to a few cluster configurations. 

Weakly bound states or resonances cannot be described in the closed quantum 
system (CQS) formalism [12]. For bound states, there appears a virtual scat- 
tering into the continuum phase space involving intermediate scattering states. 
Continuum coupling of this kind affects also the effective nucleon-nucleon in- 
teraction. For unbound states, the continuum structure appears explicitly in 
the properties of those states. 

The importance of the particle continuum was discussed in the early days of 
the multiconfigurational Shell Model (SM). However, a unified description of 
nuclear structure and nuclear reaction aspects became possible in realistic sit- 
uations only recently in the framework of the Shell Model Embedded in the 
Continuum (SMEC) [13,14] which considers the nucleus as a quantum many- 
body system coupled to the environment of decay channels (the open quantum 
system (OQS) formalism). In the SMEC, all couplings involving discrete states 
and scattering states are calculated using the SM interaction. Inclusion of the 
coupling between discrete and scattering states leads to a non-hermiticity of a 
many-body Hamiltonian that consists of the closed system with discrete eigen- 
states (the standard SM Hamiltonian), and the coupling between the closed 
system Hamiltonian and its environment of one-nucleon [13] and two-nucleon 
[15,16] decay channels. Above the particle-emission threshold, the eigenval- 
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ues of the OQS are complex. Different continuum shell-model approaches 
[17,18,19,20], including the SMEC [13,14,16], are formulated in the Hilbert 
space, i.e., they are based on the completeness of a single-particle (s.p.) basis 
consisting of bound orbits and a real continuum. A different approach to the 
treatment of particle continuum is proposed in the Gamow Shell Model [21], 
which is the multi-configurational SM with a s.p. basis given by the Berggren 
ensemble [22] consisting of Gamow (or resonant) states and the complex non- 
resonant continuum of scattering states. Gamow states and the GSM can be 
formulated rigorously in the rigged Hilbert space (see Ref. [23] and references 
quoted therein). 

In this work we present a study of the ^'^F(p,7)^^Ne radiative capture reaction 
using the formalism of SMEC [13,14]. This formalism has been applied before 
to describe the structure of mirror nuclei ^B,^Li, and capture cross sections for 
mirror reactions '^Be(p,7)'^B, '^Li(n,7)^Li [13,24] as also to the description of 
^^0(p,7)^'^F and ^^0(p,p)^^0 reactions [25,26]. The paper is organized in the 
following way. In section 2, we present a brief description of the SMEC for- 
malism emphasising those aspects of the SMEC which imply the modification 
of the CQS eigenvalues and eigenfunctions and attempts to put the theory of 
nuclear structure and reactions on the same footing. We also indicate and give 
the relevant formulae to calculate radiative capture cross sections from both 
the ground state (g.s.) and excited states of the target nucleus, with properly 
antisymmetrized wave functions being used in both initial and final states. 
Section 3 contains the results and discussions of SMEC calculations of the 
radiative capture reaction ^'^F(p,7)^^Ne. The self-consistent potential, enter- 
ing the coupled-channel (CC) equations of SMEC , and the parameters of the 
density-dependent residual interaction are presented in section 3.1. The spec- 
trum of ^^Ne obtained both in the SM and the SMEC is in section 3.2. The 
next two subsections 3.3 and 3.4 contain the astrophysical S'-factors and the 
reaction rates of the ^^F(p,7)^^Ne reaction, respectively. Finally the summary 
and conclusions are presented in section 4. 



2 Formalism 



2.1 Shell Model Embedded in the Continuum 



A unified description of interdependent nuclear structure and nuclear reactions 
is attempted in SMEC by using a projection operator technique [27]. For the 
purpose of describing the one-nucleon capture reactions, we assume that the 
Hilbert space is decomposed into two subspaces Q and P: 
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P + Q = ld 



(1) 



The Q subspace consists of A particles in many-body localized states which are 
built up by the bound state s.p. wave functions and the part of s.p. resonance 
wave functions localized inside the nucleus, i.e. before a cutoff radius Rcut- The 
P subspace contains {A — l)-particle states in localized states built up from 
(quasi-)bound single particle (qbsp) orbits and one nucleon in the scattering 
state. The rest of the s.p. resonant wave function (the Q subspace contains 
the part r < R^ut of the resonance wave function) outside the cutoff radius 
Rcut are included in the P subspace. The wave functions in Q and P are then 
renormalized in order to ensure their orthogonality in both subspaces. The 
fully antisymmetrized wave function can then be used to study the radiative 
capture processes where the capture can occur from both the g.s. and the 
excited state of the target nucleus A — 1. This method is also fully symmetric 
in treating the resonant and non-resonant part of the reaction and one does 
not need to calculate them separately. 

The first step in SMEC involves the generation of localized states in the Q 
subspace without any coupling to the continuum. This involves solving the 
standard eigenvalue problem: 



where Hqq = QHQ is identified with the SM Hamiltonian and $i are the 
localized many-body wave functions with eigenenergies Ei. 

The second step consists of solving the CC equations: 



where Hpp = PHP and stands for the outgoing boundary conditions. 
The index c = {It,l,s,j) denotes different channels characterized by It, the 
spin of the {A — 1) system, s, I and j which are the spin, orbital angular 
momentum and the total angular momentum, respectively, of the particle in 
the continuum. In a more explicit manner, eq. (3) is written as [14]: 



(2) 



- Hpp)it^ = 



(3) 







(4) 



c 



where 



(5) 
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In the above equation, Tkin stands for the kinetic-energy operator and U is the 
finite-depth potential of Woods-Saxon (WS) type with the spin-orbit term: 



U{r) 
where 



Vofir) + Vsoi^l-s: 



I dfjr) 
r dr 



c 



(6) 



f{r) 



exp 



-, -1 



(7) 



with the Coulomb potential Vc calculated for a uniformly charged sphere of 
radius Rq. v^^i in (5) is the channel-channel coupling generated by the residual 
interaction between the Q and P subspaces and depends on the target states 
of the (A — l) system, s.p. orbitals l,j and total angular momentum and parity 
J'^ of the A-particle system. The explicit formulae for H^^i and v^^, is given in 
Ref. [13]. As a residual coupling between the Q and P subspaces, we take the 
density-dependent contact force [28]: 



■■ex 
^00 



HpQ = u{p{ry^, + [l-p{rM 

+ n.r2[p(rKi + [l-p(r)]t.-]}5(ri 



1-2 



(8) 



where u is the relative strength parameter and p{r) is the nuclear density 
having a Fermi form: 



p{r) 



exp 



d 



-1 



(9) 



with r„ and d defining the size and surface thickness of the nucleus, respec- 
tively. The values taken for the parameters are mentioned in the next section. 

The diagonal part of v^^, changes the average potential U and consequently the 
s.p. wave functions. So a self-consistent procedure is adopted to calculate the 
wave functions of those s.p. states which are occupied in the target nucleus. 
This iterative procedure yields a new self-consistent average potential: 



which is used to solve eq. (4). 
The third system of equations: 



(10) 



'Ill 
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defines the functions uji tliat describe tlie continuation of the SM states $j 
into the continuum. The source term Wi couples the wave function of the Q 
subspace with those of the P subspace. 

Eigenstates of the CQS coupled to the external environment of decay channels 
are found by solving the eigenvalue problem for the effective Hamiltonian: 



H^Q^^iE) = Hqq + HqpG'^\E)Hpq (12) 

in the function space (Q) of the discrete states. Gp^\E) is the Green function 
for the motion of a s.p. in the P subspace and is given by: 

Gp = P{E - HppY^P . 

The SM Hamiltonian Hqq is hermitian. On the contrary, Htj^ is an energy de- 
pendent, no n- hermitian (complex and symmetric matrix) operator above the 
particle emission threshold and a hermitian (real) operator below the emis- 
sion threshold. The Q — P coupling contained in Hqq introduce an external 
mixing of SM states (eigenstates of Hqq) via the coupling to the decay chan- 
nels. Diagonalization of Hqq by the orthogonal but in general non-unitary 
transformation: 



= E,%*^ (13) 

yields complex eigenvalues Ei — ^iTi , which depend on the energy E of the 
particle in the continuum. The coefficients bji in (13) form a complex matrix 
of eigenvectors in the SM basis satisfying 



'^bjkbik = Sji (14) 

k 

The eigenvalues of Hqq{E) at energies Ei{E) = E, determine the energies 
and widths, Fj, of resonance states. 

The total wave function, \E^^'^°'', for the many-body problem is now obtained 
in terms of $j , , Ui and the eigenvalues of Hqq as: 



i E - Ei + {t/2)Ti 

with the incoming wave only in channel Cq as indicated by the superscript, 
and 
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(16) 



with u)j defined by Gp HpQ^i. It is to be noted here that the SMEC formahsm 
is fully symmetric in treating the continuum and bound state parts of the 
solution: v]/^^") represents the continuum state modified by the discrete states 
and Cli represents the discrete state modified by coupling to the embedding 
continuum states. 

The asymptotic conditions of Eq. (15) have been analyzed in Ref. [17]. One 
can then calculate the amplitude of the partial decay width as: 



7r = V2^ (^l^j E E / rf<^^^VX(r) (17) 

Using the proportionality relation between the matrix elements in Eq. (15) and 
the amplitudes of the partial widths (eq. (17)), one can derive the scattering 
matrix (S'-matrix): 



E-Ej + (i/2)r,- 

In the above equation S^, ^ is determined from the asymptotic, large distance 



behaviour of The amplitudes of the partial widths 7j as well as the 

complex eigenvalues Ej — {i/2)Tj of the effective Hamiltonian Hqq which 
enter in the S'-matrix, are energy dependent. 7| are also complex due to the 
channel-channel coupling and the external mixing of SM eigenstates via the 
coupling to the continuum. The partial width for channel c can then be defined 
by =1 7^ p, though the total width is no longer a sum of partial widths 
[29]: 



El7;r=f,El&.d'>f, (19) 



This is because the complex eigenvectors bji are normalized in the sense of eq. 
(14), which implies: 



E I r> 1 (20) 



and which again is a direct consequence of the non-hermitian nature of the 
effective Hamiltonian Hqq in Q. 
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Thus the formahsm of SMEC provides a way to describe nuclear structure 
and reactions in a single framework starting from the same Hamiltonian. In 
the next sub-section we shall use the total wave function (15) to describe the 
radial part of the initial state wave function in a radiative capture process. 
More details about SMEC with one particle in the continuum, and general 
discussions about the OQS could be found in Ref. [14]. 



2.2 Radiative capture 



We now give the general formulae to calculate radiative capture cross sections 
from both the g.s. and excited states of the target nucleus, with antisym- 
metrized wave functions used in both initial and final states. We shall apply 
this to the reaction ^^F(p,7)^^Ne, where the proton capture can occur from 
both If = 5/2+ g.s. and It = 1/2+ first excited state of the target nucleus, 
which is ^''F here. 

The initial wave function of the system [^''F + is: 



la,ja,It 



la ^C,ia,h^^) 



k„r 



{Y 



la 



X 



X 



(21) 



where "^i^^ describes the radial part of the wave function (15) with one 

particle in the continuum. The total angular momentum in the initial channel, 
Jj, is the result of coupling the target spin If with the total angular momentum 
ja being carried by the projectile which itself is obtained by coupling the 
intrinsic spin s of the projectile and its relative orbital angular momentum la- 



The final wave function j for the 
angular momentum, Jf, is: 



ri8 



NeY-f coupled to the final state total 



X 



s\3b 



Jf 



(22) 



where A 



hsjb 



is the coefficient of fractional parentage and u 



f 

lb,ib,hi, 



ir) is the 



s.p. wave in the many-body state Jf. s, lb and jb represent the spin, the orbital 
angular momentum and the total angular momentum of the captured nucleon 
(projectile), respectively and If^ is the final state target spin. The explicit 
summations over the target spin state in Eqs. (21) and (22) ensures that we 
take into consideration the capture from the g.s. as well as the excited states 
of the target nuclei. 
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With the wave functions ^i(r) and we can calculate the electric tran- 

sition amplitudes of multipolarity A as: 



T^^ = CiEXy- Jfhjbja{\OJjm J I J^m.;) (/,0 AO | /,0) 



X 



and the magnetic dipole transition as: 

T^'^=l'-^^jJ{l^Jfmf\Jimi){W{JatlJ^] JfJaYjJb 

+9t{-iy^-''lJtW{ItjbUf,JfIt)S,^^^^^^ (24) 



In the above formulae, 9 = rrii — rrif, a = J {2a + 1), a = Ja{2a + 1), and 



\m-^ J 

/i is the reduced mass of the system, rua and rrit are the masses, Za and Zt are 
the charges and Qa and gt are the gyromagnetic ratios of the projectile and 
the target, respectively, /ia? is the nuclear magneton and Ii'^j^j^j^, the overlap 
integral is given by: 



The Kronecker symbol 5/t7t^ in eqs. (23) and (24), ensures that the double 
summation over the target spins arising from eqs. (21) and (22) is reduced to 
a single one (over If, say). The radiative capture cross-section is then expressed 
as: 




where and ka are the wavevectors of the emitted photon and the incoming 
projectile, respectively. These formulae take into account radiative capture 
from both the excited state and the g.s. of the target nucleus. 
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3 Results and discussion 

In this section we sliall present tlie SMEC results for tlie spectrum of ^^Ne 
and tlien go on to discuss tlie astropliysical S-factor and the reaction rate for 
the ^'^F(p,7)^^Ne radiative proton capture reaction. 

3.1 The self- consistent potential 

The Q subspace in SMEC is constructed by the self-consistent, iterative method 
which for a given initial average s.p. potential (6) and for a given residual two- 
body interaction between Q and P subspaces (8) yields the self-consistent s.p. 
potential depending on the s.p. wave function /j, the total angular momen- 
tum J of the y4-nucleon system [13,14]. The channels in P: {J'^_i,lj)'^ , are 
defined by the states J^-i of the A — 1 system and the s.p. wave function Ij; 
both are coupled to J in the A-system. The parameters for the initial poten- 
tial f/(r) used in the calculations for the self-consistent potential U^'^{r) are 
Vo = 56.8853 MeV, Vso = 4.9014 MeV, Rq = 3.276 fm and a = 0.58 fm. 

1 /3 

The potential radius Rq has been defined using the convention Rq = tq^A^ ), 
where rg = 1.27 fm and At = A — 1 is the mass number of the target. In the 
self- consistent potential, we obtain the proton s.p. orbit lsi/2 at the exper- 
imental binding energy (-2.034 MeV) of the J'" = 2f level of ^^Ne and the 
s.p. orbit at the experimental binding energy (-3.92 MeV) of the g.s. 
= 0+ of i^Ne. This also puts the Opi/2 and O/7/2 s.p. orbits at -15.3 MeV 
and 8.52 MeV, respectively. Another set of parameters which enter into the 
calculations are those of the density-dependent residual interaction (8). These 
are determined to be: Vqq = 53.6,-000 = — 438.1, -Dqi = 160.9, Wof = 169.9 all in 
units of of MeV-fm^, u = 1.374, r„ = 2.64 fm and d = 0.58 fm, such that the 
experimental width [8] of the many-body state = in ^^Ne is reproduced. 

3.2 Spectrum of ^^Ne 

Fig.l presents the energy spectrum of ^^Ne calculated both in SM and in 
SMEC using the Zuker-Buck-McGrory (ZBM) interaction [30] in the Q sub- 
space. SM results with the original ZBM force [30] are shown on the l.h.s. of 
Fig.l and are labelled by 'ZBM(1969)'. With certain changes in the matrix 
elements of the effective interaction (cf Table 1) we generate the SM spectrum 
labelled by 'ZBM'. This interaction is used also to generate the SMEC spec- 
trum (labelled 'SMEC'). The experimental energy levels (labelled 'EXP') are 
plotted on the r.h.s. of Fig.l. The experimental position of the = 3f level 
is taken to be the reference level in the spectra. The SM energies , SMEC 
energy and widths, and also the experimental data are presented in Table 2. 
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Fig. 1. The spectrum of ""^^Ne. For more details, see the discussion in the text. 
Table 1 

The antisymmetrized two-body matrix elements in ZBM (this work) which differ 
from the original ZBM interaction [30]. Here, d = and s = lsi/2- 



Interaction 


[dd dd] 


[dd dd] 


[dd dd] 


[ds ds] 


[ds ds] 




T=l, J=0 


T=l, J=2 


T=l, J=4 


T=l, J=2 


T=l, J=3 


ZBM (this work) 


-2.01 


-1.61 


-0.48 


-0.57 


0.56 


ZBM (1969) 


-2.41 


-1.21 


-0.08 


-1.17 


1.16 
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Table 2 

SM energies and SMEC energies and widths vs. experimental data for ^^Ne. The 
experimental data are from: a = Ref. [31], b = Ref. [8], c = Ref. [5]. Energies and 
widths are given in keV. 



oxaxe 






SMEC 


Experiment 


TTT 
J 


energy 


energy 


energy 


width 


energy 


■ 1 j_1 

width 


n+ 
u 


— OoDD 


— 4DDo 


—4372 




—3921.6 ± 4.7 




z 


— ZODi 


OQ/1 K 

— Z04:0 


—2050 




—2034.3 lb 0.2" 




A + 
4 


— loZ 


— Doy 


—482 




—545.4 ± 0.2 




u 


1 ^9Q 

— iozy 


— D / D 


-395 




-345.3 ± 2.0'' 




2+ 


-981 


-340 


-119 




-305.2 ±0.6" 




1" 


-1326 


-644 


65 


« 


595 ± 5^^ 


0.1 ±0.1 


3+ 


600 


600 


600 


18 


600 ± 2* 


18 ±2 


0+ 


-319 


338 


941 


0.91 


666 ± 5'' 


1.0 ± 1.0 


2+ 


490 


1033 


1569 


40.2 


1182 ± 8^= 


50.0 ±10.0 


3- 


-759 


-235 


425 


0.002 


1229 ± 8^= 


< 20 


2- 


-167 


84 


839 


0.104 


1530 ± 8^= 


< 20 



3.2.1 Further effects of the continuum coupling 

Figs. 2 and 3 show the energy dependence of the eigenvalue of the effective 
Hamiltonian Hqq (eq. (12)) for J'^ = 3f and 22 , respectively. The purpose of 
these figures is not only to show the energy dependence of the eigenvalues of 
the effective Hamiltonian but also to point out that the continuum coupling 
affects both the bound and the resonant states in the energy spectra. The top 
half of the figure, in both cases, shows the real part (Eji) and the bottom 
half shows the imaginary part (Fr) of the eigenvalue as a function of the 
cm. energy (Ecm)- In all cases the solid line is the result which includes 
couplings both to the 5/2"'' g.s. and to the 1/2"'' excited state of ^^F, while 
the dashed line shows the result which includes couplings from the 5/2"'' g.s. 
only. It is interesting to note (see Fig. 2) that differs nearly by 400 keV, 
at the position of the J'" = 3f resonance, for the calculation with both the 
g.s. and the excited state of ^^F as against only with the g.s. This shows that 
the coupling to the excited state is essential for the description of the real 
part of the = 3f energy eigenvalue. The sensitivity of the width (F/j) 
of the J'^ = 3i state to the coupling to the excited state is rather weak at 
low excitation energies. At higher energies, this coupling strongly modifies the 
external mixing of unperturbed SM states, changing qualitatively the energy 
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-2 -1 1 2 3 4 5 6 
EcM (MeV) 

Fig. 2. The energy dependence of the eigenvalue (both real Eji and imaginary Tji 
parts) of the effective Hamiltonian Hqq for the 3^ state in ^^Ne. The sohd line 
corresponds to inclusion of couplings to both the g.s. and the first excited state in 
^^F, while the dashed line corresponds to the coupling to the g.s. of ^^F only. 

dependence of the width (see Fig. 2). 

A similar inference can be drawn from the energy dependence of the eigenvalue 
for the J"" = It state in ^^Ne (see Fig. 3). Inclusion of the excited state along 
with the g.s. of ^"^F results in a difference of about 150 keV at the position 
of the l\ state with respect to the calculation without the excited state of 
^^F. However, considering the relative energy shifts, including the coupling 
to the excited state seems to be more important for the ?)X than for the l\ 
state. Similarly as for the J'^ = 3f state, the imaginary part of the J'^ = 2^ 
energy eigenvalue saturates at higher Eqm in the absence of the coupling to 
the excited state of ^'''F. 

In the following, we shall study of the effect of the continuum coupling on the 
spectroscopic factors in the J'^ = 22 and J'^ = eigenvalues of the effective 
Hamiltonian Hqq. In SMEC, the spectroscopic factor, like the expectation 

value of any other operator O, can be calculated as: 



{O) 



(27) 
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Fig. 3. The same as in Fig. 2 but for the 2^ in ^^Ne. For more information see the 
caption of Fig. 2 and the discussion in the text. 

where $j (cf eq. (13)) is the eigenvector of Hqq. For the case of the spectro- 
scopic factor one identifies: 



d = a)\t){t\a (28) 

where \t) is the target state of the {A — l)-system and and a are creation 
and annihilation operators, respectively. 

In Fig. 4, we show the lsi/2 spectroscopic factors ^'^¥{b/2i) + lsi/2] — ^ 
^^Ne(2^) , (A; = 1, 2, 3) in SMEC (thicker lines) for the first three 2+ states in 
^^Ne and compare them with those obtained in SM (thinner lines). The SMEC 
spectroscopic factors are energy dependent, as a consequence of the energy- 
dependence of Hqq. Fig. 4 also shows that the coupling to the continuum 
has an overall effect of strongly modifying the s-state content of the 22 and 
23" states, increasing and decreasing them, respectively. On the contrary, the 
well-bound 2^ state is weakly affected by this coupling. This redistribution of 
SM spectroscopic factors is a clear indication of the strong mixing of a weakly 
bound 22 state and a resonance 2'^ due to the Q-P continuum coupling in 
the channel (5/2^, Si/2)^^. One may also notice that the external mixing of 2"*" 
SM states changes at higher energies close to the position of the 23" resonance 
[Er = 1.57 MeV). 
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Fig. 4. lsi/2 spectroscopic factors to the g.s. of ^^F: ["F(5/2J^)-Msi/2] ^ Ne(2^) , 
(/c = 1,2,3), for 2^ states in "^^Ne as calculated in SMEC (thicker lines) and SM 
(thinner lines). The dashed, solid and dot-dashed lines show it for 2^, 22 and 23" 
states, respectively. The change of the spectroscopic factors from the SM to SMEC 
values are indicated by arrows in each case at the experimental positions of these 
states (cf Table 2). 

Oc/5/2 spectroscopic factors [i^F(5/2^) + Oc/5/2] ^ ^^Ne(2+) , {k = 1,2) in 
SMEC and SM for the first two 2+ states in ^^Ne are shown in Fig. 5. The 
change of the spectroscopic factors by the external mixing is here less impor- 
tant and involves mainly 2f and 2^ states. For the 2^ state, the spectroscopic 
factors in SM and SMEC are both very similar and small. 



0^5/2 spectroscopic factors [i^F(l/2|) + 0^5/2] ^^Ne(2+) , {k = 1,2,3) to 
the excited 1/2^ state in ^''F are shown in Fig. 6. The channel (1/2^,^5/2)^^ 
dominates here the external mixing. The energy dependence of SMEC spec- 
troscopic factors is similar as found for Is 1/2 spectroscopic factors to the g.s. 
of ^^F (cf Fig. 4 ). Again, the redistribution of the SM spectroscopic factors 
among SMEC states concerns mainly a weakly bound 2^ state and a 2^ 
resonance. 



The continuum coupling for 22 state, as probed by the energy dependence of 
its eigenvalue and the redistribution of the spectroscopic factors, goes essen- 
tially through the coupling to the channels: (5/2j'~, Si/2)^^ and {l/2f,d5/2y^- 
The i = 2 ((i-state) coupling to the excited 1/2^ state of ^^F , changes in 
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Fig. 5. Odi/2 spectroscopic factors to the g.s. of ^^F: pT(5/2J')-h0ci5/2] Ne(2^) , 
{k = 1,2), for states in ^^Ne as calculated in SMEC (thicker lines) and SM 
(thinner lines) . For more details, see the caption of Fig. 4 and the discussion in the 
text. 

particular the weakly bound 2^ state (cf Fig. 3). This part of the continuum 
coupling for 2"*" states will be strongly suppressed in the one-proton capture 
cross-section which is sensitive mainly to low-£ {i = 0, 1) values. The role of 
the l/2i state of ^'''F for the rate of the reaction ^^F(p,7)^^Ne will be discussed 
in sect. 3.4- 

Fig. 7 shows the lsi/2 spectroscopic factor p'''F(5/2^) + lsi/2] ^^Ne(3^) , as 
calculated in SMEC and SM. The redistribution of the s-state spectroscopic 
factors among different 3+ SM states by the continuum coupling is here small 
and almost independent of the total energy of the system. 



3.3 The astrophysical S - factor for the ^'^F(p, 7/'^Ne reaction 



We now present the astrophysical S-factor for the capture reaction from both 
the 5/25'' g.s. and the l/2i excited state of ^''F. In the SMEC calculations, all 
relevant transitions from the initial {J[ = 1~, 2~, 3~, 2"*", 3"'") correlated many- 
body continuum states in ^^Ne to the final bound states (JJ = 05^2)21^2)41") 
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Fig. 6. spectroscopic factors to the excited 1/2^ state in ^"^F: 



pF(l/2|) + 0^5/2] ^ ^^Ne(2+) , {k = 1,2,3) for 2+ states in ^^Ne. The dashed, 
soHd and dot-dashed hnes show it for 2^, 22 and 2'^ states, respectively. SMEC 
values are plotted with thicker lines. For other details, see the caption of Fig. 4 and 
the discussion in the text. 



appropriate for El, Ml and E2 transitions have been included^. Since the 
position of the negative parity states in our calculations are somewhat lower 
than those experimentally observed, we have shifted the position of Jf = 1~ 
to its experimental value and have used the same shift for all calculations 
involving the negative parity states. 

In Fig. 8, we present the total S-factor (solid line) for the ^^F(p,7)^^Ne reaction 
and the separate contributions: S^i for El transitions (dotted hne), S^/i for 
Ml transitions (dashed line) and for E2 transitions (dot-dashed line). The 
8^2 in Fig. 8, which is orders of magnitude smaller than the S^;! and Sa/i, is 
multiplied by 1000. Thus, S^i and Smi components determine the reaction rate 
for the proton capture by ^'^F. At low excitation energies, Smi >'^ei and above 
the resonance the S^;! dominates. The value of the S-factor at zero energy, 
5(0), is 2130.29 eV-b and the slope, dS / dEcM\EcM=o. is -6.68 x 10^^ b. 
For Ml transitions, we find Smi(0)=1441.03 eV-b and d'&Mi/ dEcM\EcM=o—' 
6.84x10"^ b, whereas for El transitions we have S£;i(0)=689.26 eV-b and 



The subscripts i and / in and JJ, respectively, denote the initial and final 
states of the transition process. 
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Fig. 7. lsi/2 spectroscopic factor to the g.s. of ^T: ["F(5/2J^) + lsi/2] ^ ^^Ne(3n , 
for 3^ state in ^'^Ne as calculated in SMEC (thicker line) and SM (thinner line) . 
The change of the spectroscopic factors from the SM to SMEC values is indicated 
by an arrow at the experimental position of this state (cf Table 2). 



aSsi/9EcAf UeM=o=+2. 13 X 10-^b. 



We will now study the important contributions to the S-factor coming from El 
and Ml transitions. The contributions to the total S-factor due to transitions 
to the different final states of ^^Ne are shown in Fig. 9. Transitions to the 
= 2i bound state of ^^Ne dominate over all other contributions. This fact 
is also borne out by calculations presented in Ref. [10]. 

Fig. 10 shows the S£;i-factor and the individual contributions relevant for the 
capture reaction. Since the experimental Jf = 35^,25^ states in ^^Ne are at 
energies above 1 MeV, we show only their non-resonant contribution to the 
S^i-factor. The non- resonant contribution from the Jf = 3~, 2~, to JJ = 2^2 
final states are shown by thick dashed and dot-dashed lines, respectively. The 
contributions from the Jf = 1~ initial state to = 2^2 ^^^d 05^2 ^"^^ shown 
by long dashed and dotted lines, respectively. The total '^ei is represented by 
the solid line. In all CclSGS, clS has been pointed out earlier, the contributions 
to the JJ = 2i bound state of ^^Ne dominate over all other contributions. 
The low energy direct part of the El cross-section in our calculations is less 
than what has been reported in Refs. [4,10] due to a different choice of the 
potential radius. The physical picture in SMEC consists of generating the s.p. 
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Fig. 8. S-factor (solid line) for the "^^F(p,7)^^Ne reaction and the total contributions 
from the El (dotted line), Ml (dashed line) and E2 (dot-dashed line). 

wave functions of the nucleon (proton) moving in the potential well created by 
the nucleons in the target nucleus (i.e., Rq ~ 17^^^) and not in the combined 
well of the target nucleus and the nucleon (i.e., Rq ~ (17^/'^ + 1^^'^)) as used in 
[9,4,10]. This physical picture has been followed in all previous SMEC studies 
of the radiative capture reactions [13,25,26], giving an excellent description of 
the El component. 

Attempts to extract the direct capture cross section of the ^'''F(p,7)^^Ne reac- 
tion from the transfer reactions measurements are underway [32] . This proce- 
dure involves measuring the ^^N(^^F,^®Ne)^^C proton-transfer reaction cross 
section and determining the asymptotic normalization coefficients (ANC) for 
the ^^Ne bound states from a DWBA analysis of the results. Asuming that the 
measured transfer reaction is peripheral, one could then relate the extracted 
ANC with the direct capture rate at low energies [33,34]. Earlier applications 
of this method have been to the ^^0(p,7)^^F [34] and ^Be(p,7)^B [33] reactions 
to determine the astrophysical S-factor at low energies. Peripheral transfer re- 
actions would be most sensitive to the tail of the wave functions beyond the 
range of the nuclear potential than to the inner region and hence the strength 
of this method is supposed to lie in the fact that the potential model bound 
state wave function could be well reproduced by a linearly scaled Whittaker 
function in the asymptotic region [35], thereby reducing the dependence on 
the potential parameters. However it is to be mentioned here that the low 
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Fig. 9. The S-factor in SMEC is shown separately for transitions to different final 
states of ^^Ne. 

energy part of the S-factor may not be entirely due to the pure El component 
for a particular radiative capture process. Importance of this fact will become 
more apparent as we continue our discussions to the Ml transtions, especially 
on the J[ = 2+ to JJ = 2^^ Ml transitions. 

In Fig. 11, we present the contributions to Smi from J[ = 3"*" (continuum 
states) to JJ = 2^2) 45*" (bound states) transitions in ^^Ne. As remarked ear- 
lier, the transitions to the 2f final state, shown by the solid line in Fig. 11, 
dominate peaking at the position of the Sf resonance. The transitions to 22 
and to if are shown by dotted and dashed lines, respectively. The proton cap- 
ture from the excited state of ^^F has a negligible effect on the value of Smi- 
This, at first instance, may sound contradictory to the inferences from Fig. 2, 
where inclusion of the excited state couplings were shown to affect the energy 
eigenvalues. However, these couplings involve a d-state whose contribution to 
the proton capture is strongly reduced in comparison with the s-state contri- 
bution. The latter one is present in the radiative capture in the 5/2"*" (g.s.) 
channel. Moreover, the 1/2"*" excited channel opens at 495 keV, reducing the 
available energy for the proton capture and, hence, reducing importance of 
this channel in the proton capture cross-section, in particular at low Ecu- 

Much more important contribution to Sj\/i comes from transitions from the 
= 2+ continuum to the = 2'1^2 bound states of ^^Ne (see Fig. 12). The 
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Fig. 10. El part of the S-factor and the individual contributions relevant for the 
^'^F(p,7)"^^Ne capture reaction. 
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Fig. 11. Contributions to Sj\/i from J[ = 3"*" to JJ = 2f2j^t transitions in -^^Ne. 
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Fig. 12. The Ml transitions from J[ = 2+ to JJ = 2{.^ bound states of ^^Ne. 



solid line and the dashed line show the transition to the 2^ and 2^ bound 
states of ^^Ne, respectively. The 2+ continuum at low excitation energies is 
strongly correlated by the proximity of a weakly bound 22 state which in- 
duces a resonant-like feature ('resonant-halo') in the 2^^ continuum owing to 
its large s-spectroscopic factor (cf a discussion of Fig. 4). Hence, the dominant 
Ml component at low energies (below 0.4 MeV) comes from the transitions be- 
tween 22'-parented continuum states and the well-bound 2i state. The weakly 
bound 22 state plays in this scenario a role of the catalyzer of the proton cap- 
ture reaction. At higher energies, especially in the vicinity of the S"*" resonance 
at 600 keV we see the dominance of Jf = to J 



f 



2'1 Ml transitions. 



For a given strength of the continuum coupling [Q — P coupling) and fixed 
radial s.p. wave functions involved in the calculation of the source term [13,14] 
(cf eq. (11)), the rate of the proton capture in Ml transitions: = 2"^ — > 
JJ = 2f, depends solely on the position of the 22 state with respect to the 
one-proton emission threshold. In Fig. 13, we show the S^/i-factor at a fixed 
cm. energy Ecm = 50 keV for the J[ = 2^ — > JJ = 2i transitions, as a 



function of the energy (E2+) of the 22 state with respect to the one-proton 
emission threshold. The WS potential generating the s.p. wave functions is 
kept fixed in these calculations (see sect. 3.1). One can see from Fig. 13 that 
changing the energy E2+ of the 22 state from -400 keV to (the one-proton 
emission threshold) changes the Sjui-factor by two orders of magnitude. 
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Fig. 13. The S-factor for the J[ = 2+ to JJ = 2f Ml transition at a fixed energy 
EcM = 50 keV, as a function of the position (E2+) of the weakly bound 22 state in 
^^Ne with respect to the proton emission threshold. 

This astonishing effect is not due to an increased radial size of the s.p. or- 
bits, in particular of the Is-state, involved in the 2^ state since the average 
potential is kept fixed for all values of E2+ . The effect seen in Fig. 13 is a conse- 
quence of the strong interference between 22 and 23 eigenvalues of the effective 
Hamiltonian which shifts a significant part of the 23" resonance strength into 
the region of low energy continuum. This is a genuine interference effect in 
the OQS: states having a similar structure and coupled strongly to the envi- 
ronment of decay channels change in such a way that one state close to the 
emission threshold aligns its wave function with the decay channel, i.e. its 
wave function becomes more similar to the channel wave function [14,36,37]. 
A similar segregation effect has been found in Refs. [36,37] in the distribution 
of widths of the resonance states. Mixing of resonances via the coupling to the 
decay channels leads to the so-called resonance trapping, i.e. a major part of 
the decay width becomes concentrated in one (or more) resonance(s) in the 
case of one (or more) open channel(s) if the strength of the external coupling is 
sufficiently strong, whereas remaining states have negligible width and, hence, 
are trapped. In the considered example of 2+ discrete states of ^^Ne, the lsi/2 
spectroscopic amplitudes in SM states 2^ and 2^ are large, indicating a strong 
mutual coupling (strong external mixing) of these two many-body states via 
the £ = one-particle continuum. This kind of coupling may lead to the sit- 
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Table 3 

The lsi/2 spectroscopic amplitudes in the 22 state of -"^^Ne for different SM effective 
interactions. 



ZBM (this work) 


WBT [38] 


PSDFP [15] 


0.701 


0.764 


0.762 



uation that one of the 2+ eigenstates of the effective Hamiltonian closest to 
the one-proton emission threshold in the studied case) increases its s-state 
content and aligns with the decay channel [{A — 1) (g) p] = [5/2"*" ® si/2] at 
the expense of other 2"*" eigenstates. The proton capture in the s-wave then 
becomes enhanced because the correlations in low-energy continuum states 
generated by an aligned subthreshold state become enhanced by a downwards 
shift of the s-state strength from the region of 2^ resonance. This message was 
also conveyed by the conclusions drawn from Fig. 4 which shows that the cou- 
pling to the continuum has the overall effect of strongly modifying the s-state 
content of the 2^ and 2^ states, and thereby redistributing their strengths. 

The effect is all the more accentuated by the large s-state spectroscopic am- 
plitude in the 22 state. To verify that this is not an artifact of the ZBM 
interaction, we have also shown in Table 3 the corresponding numbers for two 
different SM effective interactions, WBT [38] and PSDFP [15], also used for 
A = 18 mass region. We see that the magnitudes of the s-state spectroscopic 
amplitudes are similar for all these interactions. 



3.4 The -'^'^F(p, 7)^^Ne reaction rates 

As has been mentioned earlier the radiative proton capture by ^''F plays an 
important role in the nucleosynthesis. In Table 4, we present the reaction 
rate, Na{(tv) in cm^mole~^s~^ for this reaction as a function of some typical 
stellar temperatures T (in GK). (av) is the reaction rate per particle pair [39] 
and Na is Avogadro's number. Column 1 of Table 4 shows the SMEC results 
calculated taking into consideration the capture from both the 5/2"*" g.s. and 
the 1/2+ excited state of ^^F, while the column designated by SMEC (g.s.) 
shows the result for proton capture from the 5/2+ g.s. of ^^F alone. We do not 
see any noticeable difference in the reaction rates in these two calculations, 
particularly below temperatures of 0.4 GK which are more important for nova 
nucleosynthesis. This is understandable because the Gamow window spans 
lower energies at these temperatures and, moreover, the 1/2+ excited state 
channel involves the proton capture in the d-wave and is strongly disfavored 
by the centrifugal barrier. 
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Table 4 

^''F(p,7)^^Ne reaction rates in cm'^mole~"'^s^^ in SMEC and Ref.[8] for some typical 
stellar temperatures (in GK). For other descriptions see the text. 



T(109K) 


SMEC 


SMEC (g.s.) 


Ref. [8] 


0.1 


1.317 


xlO-^ 


-9 


1.317 


xlO-^ 


-9 


(2.68±0.38) xlO-9 


0.2 


2.209 


xlO-^ 


-6 


2.209 


xlO-^ 


-6 


(5.15±0.75) xlO-6 


0.3 


7.736 


xlO" 


-5 


7.736 


xl0~ 


'5 


(1.97±0.29) xlO"^ 


0.4 


7.458 


xlO" 


-4 


7.453 


xlO" 


-4 


(2.29±0.40) xl0~3 


0.5 


4.124 


xlO-^ 


-3 


4.116 


xlO-^ 


-3 


(1.77±0.49) xl0~2 


0.6 


1.646 


xlO' 


^2 


1.640 


XlO" 


-2 


(9.29±3.28) xlO"^ 


0.7 


5.023 


xlO-^ 


~2 


5.000 


xlO-^ 


-2 


(3.32±1.30) xlO-i 


0.8 


1.227 


xlO" 


1 


1.221 


XlO" 


-1 


(8.80±3.61) xlO-i 


0.9 


2.516 


xlO' 


^1 


2.503 


XlO" 


-1 


(1.88±0.78) xlO° 


1.0 


4.516 


xlO' 


^1 


4.491 


XlO" 


-1 


(3.43±1.44) xlO° 


1.5 


2.643 xlO 





2.629 


XlO 





(1.97±0.78) xlO^ 


2.0 


6.185 


xlO 





6.155 


XlO 





(4.62±1.64) xlO^ 



The last column in Table 4 is the reaction rate from Ref. [8], which is quoted 
here for the purpose of comparison with our reaction rate. The SMEC rates are 
lower by a factor of 2-3 below 0.4 GK and by a factor of 7.5 at 2.0 GK, as com- 
pared with those of Ref. [8], primarily because of the lower El contributions 
due to the reasons mentioned in the previous section. 

In Fig. 14, we show the variation of the total reaction rate as a function of 
temperature and for some important transitions. The total contribution is 
shown by the solid line, while the dotted and short-dashed lines represent the 
contributions from the El and Ml components, respectively. Amongst the Ml 
components the contributions from the J[ = 2"*" to JJ = 2^2 components 
(long-dashed) dominate the rate at lower temperatures and accounts for al- 
most all the Ml cross-section below about 0.4 GK, while the contributions 
from J[ = to JJ = 2^2 (dot-dashed) are dominant at higher temperatures. 

It would be interesting to investigate the astrophysical implications of our 
^^F(p,7)^^Ne reaction rate, especially in nova explosions. The previous reaction 
rates of Wiescher et al. [3] are higher than that of Bardayan et al. [8] based on 
the recent determination of the energy and the total width of the J'^ = state 
in ^^Ne and adding the non-resonant direct capture rate from the potential 
model in Ref. [4]. The reaction rates of Sherr and Fortune [40] and Garcia 
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Fig. 14. The variation of the total reaction rate as a function of temperature and for 
some important transitions. The total reaction cross-section is shown by the solid 
line, while the dotted and short-dashed lines represent the El and Ml components, 
respectively. For more details, see the description in the text. 

et al. [4] are rather similar to those of Bardayan et al. below about 0.4 GK 
and are a bit lower at higher temperatures. The reaction rate of Ref. [8] was 
used to calculate [41] the nucleosynthesis in nova outbursts on the surfaces of 
1.25 Mq and 1.35 Mq ONeMg white dwarfs and a 1.00 Mq CO white dwarf. 
It was shown that for a 1.25 Mq white dwarf a slower -'^'^F(p,7)^^Ne reaction 
rate produced more ^^F than the faster rate because the ^^F (produced from 
the /3-decay of ^^Ne) at lower temperatures would survive as ^^F and would 
not be destroyed by the ^^F(p,a)^^0 reaction. Although a faster ^''F(p,7)^^Ne 
reaction rate does allow more conversion of ^^F to ^^F, but since this ^^F would 
be produced at higher temperatures it would be destroyed by the ^®F(p,a)^^0 
reaction. Thus a larger final abundance of ^^F would be obtained with a slower 
"^'^F(p,7)^^Ne reaction rate. It was also calculated that for a 1.35 Mq white 
dwarf the abundances of ^^O and ^^F increased by a factor of 14000, as against 
estimates based on rates from Ref. [3], in the third hottest zones of the nova. 
For a 1.00 Mq CO white dwarf, which involves lower peak temperatures than 
the other two cases, only small variation in the mass fraction of various isotopes 
were observed except for ^^Ne; where it increased by 21% as compared to Ref. 
[3]. 

Our reaction rate is lower than those of Ref. [8] by nearly a factor of 2-3 
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below 0.4 GK. This is also about the temperature range of the hottest zone 
in 1.25 Mq and 1.35 Mq white dwarfs, which are 0.333 GK and 0.457 GK, 
respectively [41]. For a 1.25 Mq white dwarf, with our lower rate, we would also 
qualitatively expect a large ^^F abundance produced at lower temperatures 
and not being destroyed by the ^^F(p,a)^^0 reaction due to the same reasons 
as mentioned in the previous paragraph. The final abundances of ^^F, in 1.35 
Mq white dwarf, could be slightly less than in Refs. [3,41], following the same 
trends as in Ref. [41]. Moreover, recent measurements [42] seems to suggest 
that the estimated ^^F(p,a)^^0 reaction rate could be 2-5 times less than what 
was estimated in recent calculations [43,44,45]. This would imply the survival 
of ^^F (Ti/2 = 110 min) and this opens up the possibility of detecting the 
characteristic 511 keV 7-ray, from the electron-positron annihilation following 
the /3-decay of ^^F, using the techniques of gamma-ray astronomy. This, in 
turn would provide us an opportunity to gain further insights about the nova 
mechanism. 

Furthermore, our lower ^^F(p,7)^^Ne reaction rate would allow more ^"^F to 
survive and the production of more ^^O by ^^F(/3"'", i/)^^0 in both 1.25 Mq and 
1.35 Mq white dwarfs as compared to the values reported in [41]. Moreover, 
apart from the branch which produces ^^F via the ^^0(p,7)^^F reaction, the 
synthesis of ^^O by the ^^0(p, a)^^N(p, 7)^^© reaction chain would be more 
feasible with the rates reported in this work. The /3-decay of ^^O to ^^N could 
be a way of explaining the overabundance of ^^N in the nova ejecta [46]. 

For a 1.00 Mq CO white dwarf, which involves lower peak temperatures, we 
would generally concur with the scenario predicted in Ref. [41] as our rates 
would be closer to the rates of Ref. [8] at lower temperatures. We would 
probably have similar abundances for most of the isotopes as in Refs. [3,41] 
except for abundances of ^^Ne, which could be higher in our case following the 
similar trends as in Ref. [41]. 



4 Summary and conclusions 

In this work, we have applied the open quantum system formalism of the 
SMEC for the microscopic description of the ^®Ne spectra, low energy as- 
trophysical S-factor and reaction rates for the radiative proton capture reac- 
tion ^^F(p,7)^'^Ne. SMEC is a model in which the realistic SM solutions for 
(quasi-)bound states are coupled to the environment of one-particle [13,14] 
and two-particle decay channels [16] for the description of complicated low 
energy excitations of weakly bound nuclei and various nuclear reactions and 
decays involving up to two nucleons in the scattering continuum. This theo- 
retical model, which attempts a unified description of interdependent nuclear 
structure and reaction theories by using a projection operator technique [27], 
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is a recent development of the Continuum Shell Model (CSM) [17,18,19]. 



Inclusion of the target excited state was shown to have an influence on the 
eigenvalues of the effective Hamiltonian. It was seen that the energy of the 
resonance in ^^Ne could be different by as much as 400 keV, depending on 
whether the excited state couplings were included or not. It was also interesting 
to observe the energy dependence of the spectroscopic factors for the first three 
2+ states in ^^Ne. The coupling to the continuum strongly modifies the s-state 
content of the weakly bound 22 and resonant 23", decreasing the latter and 
reinforcing the former, with respect to their SM values. On the other hand the 
s-state content of the 2f state was not affected too much by the continuum 
couplings by virtue of it being well below the particle emission threshold. It 
would ve very nice if, in some future experiment, these spectroscopic factors 
could be measured, especially those of the 2f and 22 ■ 

We also calculated the astrophysical S-factors for the radiative capture reac- 
tion ^^F(p,7)^^Ne, where the capture could occur from both the 5/2"*" g.s. and 
the 1/2"*" first excited state of ^^F. Fully antisymmetrized wave functions were 
used in both the initial and final states of the radiative capture process. This 
method is also fully symmetric in treating the resonant and non-resonant part 
of the reaction and one does not need to calculate them separately. We have 
shown that in the region around 600 keV the J[ = 3"^ to JJ = 2f compo- 
nent of the Ml dominates the S-factor, but at very low energies (below 0.4 
MeV) the Ml component arising from the J[ = 2+ to JJ = 2f transition 
has a substantial impact on the low energy S-factor. In fact its contribution 
is even more than the El component below 0.1 MeV. This has been traced 
back to the strong mixing of the weakly bound 2^ state with the 2^ resonance 
in the low energy correlated 2~^ continuum wave function. One should, thus 
be careful while analyzing the low energy direct capture part of the reaction 
with theories, which tend to minimize the effect of nuclear structure in nuclear 
reactions. It would indeed be extremely interesting if in proposed experiments 
[47], the low energy El and Ml components could be disentangled and the 
low energy behaviour of the Ml is studied in more details. 

We have also calculated the reaction rate for the ^^F(p,7)^^Ne radiative cap- 
ture reaction for some typical stellar temperatures. Our rates are lower by a 
factor of 2-3 below 0.4 GK and by a factor of 7.5 at 2.0 GK, as compared with 
those of Ref. [8]. Qualitatively we predict a similar situation in novae nucle- 
osynthesis, following the trends reported in Ref. [41], with perhaps a larger 
abundance of ^^F and ^''O. It would also be interesting to investigate if our 
smaller reaction rates could eventually lead to a larger ^^O abundance whose 
/3-decay to ^^N could be a way of explaining the overabundance of ^^N in the 
nova ejecta [46]. 
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